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Abstract 

A theoretical description of time correlation functions for electron properties in the presence of 
a positive ion of charge number Z is given. The simplest case of an electron gas distorted by a 
single ion is considered. A semi-classical representation with a regularized electron - ion potential is 
used to obtain a linear kinetic theory that is asymptotically exact at short times. This Markovian 
approximation includes all initial (equilibrium) electron - electron and electron - ion correlations 
through renormalized pair potentials. The kinetic theory is solved in terms of single particle 
trajectories of the electron - ion potential and a dielectric function for the inhomogeneous electron 
gas. The results are illustrated by a calculation of the autocorrelation function for the electron field 
at the ion. The dependence on charge number Z is shown to be dominated by the bound states of 
the effective electron - ion potential. On this basis, a very simple practical representation of the 
trajectories is proposed and shown to be accurate over a wide range including strong electron - ion 
coupling. This simple representation is then used for a brief analysis of the dielectric function for 
the inhomogeneous electron gas. 

PACS numbers: 05.20.Dd, 45.70.Mg, 51.10.+y, 47.50. +d 



I. INTRODUCTION 



Electron dynamics in a rigid uniform neutralizing background is a well-studied problem 
(jellium in quantum mechanics-, one component plasma in classical mechanics^). More 
realistically, point ions lead to a polarization of the electron density (e.g., in a hydrogen 
plasma) and the dynamics of the non-uniformly distributed electrons is radically changed. 
The objective here is to provide a practical theory for the description of equilibrium time 
correlation functions for electrons in the simplest case of a single point ion of charge number 
Z . If the charge is positive, essential quantum diffraction effects must be accounted for even 
at high temperatures and low densities to avoid the electron - ion Coulomb singularity. A 
classical Hamiltonian description is used here, with a regularized electron - ion interaction 
that accounts for such effects^. This study is an outgrowth of recent investigations based on 
molecular dynamics simulations for this system^. The qualitative features observed for the 
electron field autocorrelation function from simulation were captured by a simple mean field 
kinetic theory. Such a kinetic theory is obtained here from the asymptotically exact short 
time limit for the generator of the dynamics, providing both context and a generalization of 
the analysis in reference 4 to strong electron - electron coupling conditions. 

The kinetic theory is solved exactly to express the correlation functions in terms of ef- 
fective single electron trajectories about the ion and collective excitations via a dielectric 
function for the non-uniform electron fluid. For Z = the results reduce to the familiar 
random phase approximation (RPA) with local field effects (the generalized Vlasov ap- 
proximation of reference 5; see also reference 6 for a related nonlinear kinetic equation for 
dusty plasmas). More generally it constitutes a generalization of the RPA to a non- uniform 
electron gas, with both ion - electron and electron - electron interactions renormalized by 
correlations (the Vlasov equation for an electron gas in a periodic potential is discussed 
in reference 7). The only required input is the time independent correlations for one or 
two electrons and the ion. For the calculations here the hypernetted chain approximation 
(HNC) integral equations are used for these static correlations. The correlation functions 
are further decomposed into contributions from the bound and free (positive and negative 
energy) states of the effective single particle dynamics in Section IIIII 

As a special case, the electric field autocorrelation function is considered in Section IIVI 
with the objective of providing a clear interpretation for the Z dependence observed in 
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simulations. For increasing Z this dependence includes 1) an increasing covariance of the 
field (initial value of the correlation function), 2) a decreasing correlation time, and 3) the 
development of a strong domain of anti-correlation at intermediate times^. It is shown here 
that all three features can be attributed to an increasing contribution from the bound states 
of the single particle effective dynamics representing actual metastable trapped trajectories 
of the particle dynamics. With this understanding of the active mechanisms, a simple 
analytic and accurate model for the bound and free state contributions is proposed and 
tested. The dynamics is restricted to circular and straight line trajectories, and the electron 
- ion charge correlation is represented by a nonlinear Debye distribution. Somewhat sur- 
prisingly, the model reproduces all of the above Z dependencies with remarkable accuracy. 
This provides the basis for a practical representation of more general correlation functions, 
such as the dynamic structure factor, and more complex state conditions required for plasma 
spectroscopy in hot, dense matter-i^. 

To illustrate the practical utility of the model, collective excitations are explored briefly in 
Section |VT] using the model to evaluate the dielectric function for this nonuniform electron 
distribution about the ion. For weakly nonuniform conditions (small Z) the results are 
suggestive of a local density approximation whereby the modes are similar to those of a 
uniform electron gas, but with the density replaced by the actual local density near the ion. 
However, this simple approximation fails for larger Z where the bound states dominate and 
long wavelength plasmons are replaced by local excitations at the circular orbit frequencies. 
Finally, some future directions are discussed in the last Section. 

II. CORRELATION FUNCTIONS AND MARKOVIAN APPROXIMATION 

Consider a system of electrons of charge — e, an inflnitely massive positive ion of 
charge Ze placed at the origin, and a rigid uniform positive background for overall charge 
neutrality contained in a large volume V . The Hamiltonian is 



where and Vq, are the position and velocity of electron a. The Coulomb interaction 
between electrons a and 7 is denoted by Veeija'^^ where r^^ = \va — Also, Vei (r^) is the 
electron-ion interaction for electron a, and V^;, (r^) is the Coulomb interaction for electron 
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a with the uniform neutrahzing background. In a quantum description, Vei [To) is also a 
Coulomb interaction but in the classical case the short range attractive divergence must be 
"regularized" within a distance 5 of the order of the de Broglie wavelength^. The simplest 
such form isi^ 

V;.M = -— (2) 

Toe 

In the remainder of this presentation such a semi-classical description is assumed. Comments 
on the corresponding quantum analysis are given in the final Discussion section. 

The typical response functions characterizing dynamical excitations in a plasma are the 
charge density or current autocorrelation functions, which are sums of single particle func- 
tions. More generally, the correlation functions of this type are defined by 

CAB{t) = {A{t)B) = j dVA{Vt)B (r) pe (r) (3) 

where F = {xi,..,Xiv^} is a point in the QN^ dimensional phase space, and Xa = rQ,,VQ, 
denotes a point in the phase space of particle a. The notation Tt denotes the evolution of 
the point F to a time t later under the dynamics generated by the Hamiltonian of ([1]). The 
role of the central fixed ion is suppressed in this notation, and it acts as an external potential 
for the electrons. The phase functions ^(F) and B{T) denote some observables of interest, 
composed of sums of single particle functions. 

A = Y^ a(a;«), ^ = h{xc?j. (4) 

Finally, the average is over an equilibrium ensemble (e.g., Gibbs), Pe (r)- Because of the 
special form the particle average can be reduced to a corresponding average in the 
single electron subspace, by partial integration over Ne — 1 electron degrees of freedom (see 
Appendix [B|) 

C'ab(^) = J dxn{r)(f) (v) a{x)b(x,t). (5) 

Here, n(r) is the equilibrium number density for electrons at a distance r from the ion (the 
precise definition as a partial integral of pe (F) is given in Appendix A), and (f){v) is the 
Maxwell-Boltzmann velocity distribution. The function b{x, t) at t = is linearly related to 
the single particle phase function b{x) in (jl]) 

b{x, 0) = b{x) = b{x) + / dx'n{r')(f) {v') h (r, r') b{x'). (6) 



The correlation function h (r, r') is related to the joint number density n{r, r') for two elec- 
trons at r and r' with the ion at the origin by 

n{r)n{r')h (r, r') = n{r, r') — n{r)n{r'). (7) 

The precise definition for ?7,(r, r') as a partial integral of pe (r) is given in Appendix |X1 The 
time evolution of t) in the single particle phase space is governed by a linear equation 
of the form 

S.KM) + /<*.'£(...';«)K.'.«)^o. (8) 

All of the results up to this point are still exact. 

The difficult many-body problem is encountered in the determination of C (x, x'\ t). Weak 
coupling and perturbation expansions are not appropriate for high Z ions or conditions for 
strongly coupled electrons so instead a Markovian approximation is proposed, 

C (x, x'; t) ^ C (x, x'; t = 0) = £ (x, x) . (9) 

This approximation assumes that the exact generator for the initial dynamics persists as 
the dominant form for later times as well. In this way the exact initial correlations among 
electrons and with the ion are included. The detailed form for C (x, x') is obtained in 
Appendix [Bl with the result 

C (x, x') = (v ■ Vr - m- VrV,e (r) ■ Vv) (5 (x - x') + V ■ Vr/9Vee (r, r') iy') n ir') , (10) 

where Vie and Vee (r, r') are "renormalized" electron - ion and electron - electron inter- 
actions 

V,e (r) = -13-^ In n (r) , Vee (r, r') = -r (r, r') . (11) 
The direct correlation function c (r, r') is defined in terms of h (r, r') by 

c (r, r') = h (r, r') — j dr"h (r, r") n [r") c (r", r') . (12) 

At Z = this becomes the usual Ornstein - Zernicke equation^. 

To interpret (|TOl) . substitute this approximation into (IHl) to get the Markovian linear 
kinetic equation for 6(x, t) 

(9t + V ■ Vr - m- VrV.e (r) • Vv) 6(x, t) = -V ■ Vr/3 j dx'Vee (r, v') {v') 71 (r') h{x\ t). 

(13) 



At weak electron - electron coupling Vee (i", i"') Kjedr — r'l) and at weak electron - ion 
coupling Vie (r) — > Vie (r) , and f|T3l) is recognized as the linear Vlasov equation. More 
generally, the Markov approximation (1131) upgrades this mean field result to include the 
effects of equilibrium correlations on all interaction potentials. Thus it is suitable for a 
discussion of the strong coupling conditions that occur for Z > 1. The left side of f ll3p 
describes single electron motion about the ion in the effective potential Vie, while the right 
side describes dynamical screening of this motion. 

In summary, the description of electron dynamical correlations and fiuctuations has been 
reduced in the Markovian approximation to 

CABit) = j dxn{r)(t){v)a{x)e-'^%{x), (14) 

where C is the operator whose kernel is ffTOj) . This operator requires as input the equilibrium 
electron density n (r) and the equilibrium direct correlation function. The kinetic equation 
can be solved exactly in terms of the single particle trajectories about the ion and dielec- 
tric function for an inhomogeneous electron gas, describing the dynamical screening due to 
interactions among the electrons in the presence of the ion. The details are carried out 
in Appendix [Cl and the correlation functions are obtained from that solution in Appendix 
Id] For the class of correlation functions for which a{x) = a(r) (i.e., is independent of the 
velocity) the Laplace transform of f|T4l) takes the simpler form 

dte-'^CABit) = J drdvn{r)(t){v)a{v;z)goiz)b{x), (15) 

The dynamics is governed by the resolvent operator 

Goiz) = {z + Coy\ £0 = V ■ Vr - m- VrVie (r) ■ Vv. (16) 

The generator for the dynamics, Co, is seen to be that for a single electron interacting with 
the ion via the effective mean field potential Vie (?")■ The function a(r; z) is the given function 
a(r), modified by dynamical screening 

(17) 

where e (r, r'; z) is the "dielectric function" for the electrons in the presence of the ion^i 

5{r-r')- j dr"7r (r, r"; z) Vee{r", r'), (18) 



and 71 (r, r"; z) is 

—j3n{r 



J tiv0(t;)6?o(^)v-Vr5(r-r"). (19) 
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A. Dynamic structure factor 



An important example is the autocorrelation function for the electron density near the 
ion. In the absence of the ion this is referred to as the dynamic structure factor and that 
terminology will be used here in the presence of the ion as well. The correlation function 
CAsit) = C(q, q';t) is constructed from the local densities of (jl]) with a{xa) = ^ (q — 
and b{xa) = 5 (q' — ra). Then f|T5|) becomes 

j rfte-"*C(q, q'; t) = J drdwn (r) {v) e"^ (q, r; z) Go (z) s (r, q') . (20) 

Here s (r, q') is the static structure factor 

s (r, q') = 5 (r - q') + n {q') h (r, q') = 5 (r - q') + {q') h (q', r) 

= e-i(q',r;0), (21) 

representing the exact initial correlations. The last equality of (^T^ is proved in Appendix 

El 

Equation fl20|) is the exact short time (Markovian) form for the dynamic structure factor. 
For Z = it becomes the usual random phase approximation (RPA) with "local field 
corrections" ; this means that the bare electron - electron potential has been replaced by Vee 
(the corresponding direct correlation function) to account for exact initial correlations. The 
Z = case has been studied in detail for the hydrogen plasma, where this approximation is 
found to be very good up to moderate plasma coupling strengths over a wide range of space 
and time scales^. Equation (fT6l) extends this approximation to include the presence of the 
ion for Z ^ 0, corresponding to an inhomogeneous RPA. 

B. Electric field autocorrelation function 

A second important example is the autocorrelation function C [t) for the electron electric 
field at the ion, where 

A = B = J2^i^a), e(r,) = Vr„Ki(r«). (22) 

a=l 

This correlation function also is obtained from the dynamic structure factor C(r,r';t) by 
integration 



so 
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^ dte~''C{t) = J drdvn (r) (v) e(r; z) ■ Go (z) e,(r). (23) 
Interestingly, one of the fields is dynamically screened while the other is statically screened, 

e(r; = j dr'e{r')e-^ (r', r; z) , e,(r) = J dr'e{r')s (r', r) = e(r; z = 0). (24) 

III. BOUND AND FREE CONTRIBUTIONS 

The trajectories of the mean field generator Cq do not represent the dynamics of any 
given electron, but rather their effective collective representation. Thus, the interactions 
among many electrons appears in the mean field theory only through their modification of 
the potentials. The bound states of the effective potential Vie (f) are representations of real 
metastable states in the MD simulation. At weak coupling there are few such metastable 
states and their lifetimes are short compared to the correlation time for the field autocorrela- 
tion function. As Z increases the stronger coupling gives rise to more metastable states with 
longer lifetimes. There is a crossover of these lifetimes to values larger than the correlation 
time at which point they behave essentially as bound states for the relevant time scales. 
To isolate the effects of such bound states it is useful to divide the phase space integral of 
f|T^ into contributions from bound and free parts of that phase space. The decomposition 
is defined by the negative and positive energy states for the effective potential Vie ('")• For a 
given position r there is a maximum velocity Vm{r) above which the total energy is positive 

vUr) = ^-2Vie (r) /m. (25) 

The single particle equilibrium density for the position and velocity can therefore be divided 
into two contributions 

n(r)0 (t>) = 6 (fm(?") — "w) n(r)0 (t>) + 6 (t; — fm('")) '^(?")0 (f ) 

= {n{r)4>{v)\\{n{r)4>{v))^. (26) 

For example, integration over the velocity gives the relative contribution of bound states to 
the density n(r) 

nb(r) = / dv {n{r)(j){v))f^ 



8 



Jo 

= n{r) ferf ( - ^ . (27) 

where Vq = \j2kBT /mf. is the thermal velocity of the electron and erf (x) is the error function. 
Since Vm{r) is proportional to \/Z, nh{r) is an increasing function of Z for all r. 

The decomposition (126!) provides the identification of contributions to the correlation 
functions from bound and free states 

CAB{t) = CUt)+CUt)- (28) 

The analysis of the following Sections shows that the interesting Z dependence of electron 
dynamics can be understood in terms of the relative sizes of these two contributions. 



IV. EXAMPLE: ELECTRIC FIELD AUTOCORRELATION FUNCTION 

The dielectric function e (r, r'; z) describes a crossover from no screening at short times 
(large z) to static screening at large times {z 0). In the next two sections, the short time 
form with e (r, r'; z) 6 {r — r') will be considered. In that case the correlation functions 
in (|T5|) become the effective single particle functions 

CAsit) ^ y rfrdvn(r)0(w)a(r)e-^°*6(x), (29) 

which are asymptotically exact at short times. The effects of dynamical screening at longer 
times are discussed briefly in Section IVll 

The electric field autocorrelation function is particularly instructive since the field is 
sensitive to configurations closest to the ion. Also, its time integral determines the domi- 
nant contribution to the half width of spectral line widths broadened by electrons in many 
practical cases^^. The short time form of fl23l) is 

C{t) j drdvn (r) (v) e(r) ■ e~^«*e,(r). (30) 

It is shown in Appendix |A] that the statically screened field of ( l24l) simplifies further to 

e,(r) = ^VrV.e(r). (31) 
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FIG. 1: Bound and free state contributions to the field covariance C(0) as a function of Z. 

Thus, the number density n (r) determines all of the ingredients needed for calculation of 
Cif). It is calculated here in the HNC approximation described in Appendix |Al The electron 
- electron coupling strength is measured by the dimensionless ratio F = /3e^/ro where tq is 
the average distance between electrons, determined from the density by ^'nnj-^j?) = 1. The 
ion - electron coupling is measured by cr = —(3Vei (0) = (3Ze^ /5 = ZT (ro/5). The results 
presented here are for F = 0.1, and cr = 0.25Z for values of Z < 40. The corresponding 
quantum regularization length is 6/rQ = 0.4. The electron - electron coupling is therefore 
weak, but the ion - electron coupling can be very strong, a < 10. These conditions were 
chosen because previous molecular dynamics studies have been performed at these values^. 

It is useful to anticipate the increasing role of bound states with increasing Z by consider- 
ing first the time independent covariance C(0). This is shown in Figure[T] The sharp increase 
above Z ~ 5 is seen to be entirely due to the appearance of the bound states. Similar strong 
effects on dynamical structure are observed. Figure [2] shows the results for C(t) calculated 
from ( l30i) for Z = 1,4,8,20,30,40. The development of a strong anti-correlation and the 
decreasing initial correlation time with increasing Z is evident. These are the effects noted 
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FIG. 2: Autocorrelation function for an electron field at an ion of charge number Z = 1, 4, 8, 20, 30, 
and 40 

above, first observed in MD simulations^. The interpretation of this two-fold dependence 
on Z is provided by Figures [3] and H] showing the contributions from bound and free state 
contributions for Z = 4, 30. For Z = 4 the dominant contribution is from free states, which 
have a monotonic positive decay. In contrast, for Z = 30 the dominant contribution is from 
bound states which provide the negative anti-correlation as the sign of the field changes 
along each trajectory when it passes through apsidal distances. The time for this change 
can be estimated by half the period for a circular orbit at position r, which is proportional 
to (r^/Z) . This is consistent with the observed decrease in correlation time in Figure El 
Further elaboration and explanation is provided by the simple model of the next section. 

The results of this section demonstrate the utility of the kinetic theory for conditions of 
strong ion - electron coupling. Although only weak electron - electron coupling was con- 
sidered, the theory is applicable to strong coupling among electrons as well. Also, while 
attention in this section has been limited to the electric field autocorrelation function it 
is clear that the analysis apphes with equal ease to the dynamic structure factor as well. 
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FIG. 3: Bound and free state contributions to C{t) for Z = 4. 

with only the additional complication of more parameters (i.e., two position vectors) char- 
acterizing that function. The decomposition of the correlation function into bound and free 
contributions demonstrates that the interesting features associated with increasing Z can 
be attributed entirely to the increasing contribution from bound states. 



V. A SIMPLE, ANALYTIC, AND ACCURATE MODEL 

Consider again the electric field auto correlation function given by fl30|) explicitly decom- 
posed into its bound and free contributions 

C{t) = C\t) + Cfit), (32) 

C^'^(t) ^ j drdw {n (r) {v))^^ e(r) ■ e-^°*e,(r), (33) 

The objective here is to capture the qualitative features of the bound and free contributions in 
a very simple model that allows further elaboration of their relative roles and the mechanisms 
involved. This is accomplished by assuming circular trajectories for the bound states and 
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FIG. 4: Bound and free state contributions to C{t) for Z = 30. 

straight line trajectories for the free states, 

{t) / drdv (n (r) {v))j e(r) ■ e,(r - vt), 



(34) 
(35) 



drdv (n (r) (f ))b e{r)es{r) cos 

Here cos(fc(r)t/r) = r ■ r{t) and r(t) = r for circular obits. In this case the velocity must 
be orthogonal to r with the specified magnitude 



r dVie (r) 



(36) 



m dr 

for consistency with Newton's equations. The Maxwellian in (n (r) (f ))^ must therefore be 
replaced by this restriction on the velocities 



[n (r) {v))f^ —Uf, (r) {vm{r) - v) 6 (v-r) 6 {u - (r)) 



(37) 



where u is the component of v orthogonal to r. The factor ni,{r) is given by fl271) and is 
required by the correct normalization on integration over all velocities. Equation (!35l) then 
becomes 

POO 

C\t)^A7r drr^nb{r)e{r)es{r)cos{uJc{r)t) (38) 
Jo 
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Kinetic theory and simple model 
Z=4,30 r=0.1 a=0.25Z 



theory (solid lines) 
model (circles) 
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FIG. 5: Comparison of C{t) calculated from kinetic theory and from the simple model for Z = 4 
and 30. 

where uJc{r) = Vc{r)/r. Note that these modifications of the trajectories do not affect the 
initial value, C(0), which is still exact. 

One further simplification is made to complete the model. The effective potential Vie {^) 
is replaced by its weak coupling Debye form 

Ze^ 1 



V.e (r) 



-r/5 



). 



(1 - (S/Xf) r 

where A = ro/y/ST is the Debye length and 6 is the quantum regularization length of 
The corresponding electron density is now the non-linear Debye form 



(39) 



n(r] 



Ze^ 



.i(e-/A_e-r/5 



(40) 



(1 - {5/XY) r 

This is exact in the weak coupling limit, with Z ^ Z. More generally, Z is chosen to give 
the correct value of C(0) using (!39l) and (HOl) in the exact equation for C(0) and adjusting 
Z to fit the values obtained from HNC 



C*(0) = ^C(O) = 4- / dyye-y (e'f^y^ - 1 
ZT Jo ^ 



(41) 
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with 

^s- (i- (vsTi-) y ' 

Table 1 gives the values obtained using C*(0) from the HNC approximation, for the case 
r = 0.1 and a = 0.25Z. 



z 


Z 


C*{0) 


1 


1.03 


2.67 


4 


3.76 


3.13 


8 


7.29 


3.92 


20 


16.79 


8.15 


30 


22.87 


14.4 


40 


27.20 


22.7 



Table 1: Effective charge number Z for the Debye form 

Figure 5 shows again the results of Figures 3 and 4 for Z = 4 and 30, now including as 
well the results from the simple model of this section. Remarkably, the use of the Debye 
form with the straight line and circular trajectories gives an accurate representation of the 
kinetic theory results. This provides the basis for a practical tool for use in more complex 
conditions, as discussed in the last section, and for other correlation functions, as illustrated 
in the next section. 

VI. DIELECTRIC FUNCTION 

The last two sections addressed the short time form ( l29i) of the correlation functions 
for which the dielectric function behaves as e (r, r'; 2;) — > 5(r — r'). For very long times 
(small z) e (r, r'; z) crosses over to represent static screening. At intermediate times there 
are contributions from collective excitations. Their description is more complex than for the 
uniform electron gas for which the space dependence of the dielectric function occurs only 
through r — r'. To simplify the discussion here consider the weak electron - electron limit 
(but possibly strong ion - electron coupling) for which Vee(r, r') — * Vee(|r — r'|) and define 
the partial transform 

15 



dte-'' J dw(j) {v) ik ■ V (-t) e^'^-('-(-*)-'-) . (43) 

Collective excitations for the non-uniform system, z (k, r), are defined by ?(r, k; z (k, r)) = 0. 
For zero charge number on the ion, Z = 0, ?(r, k; 2;) reduces to the familiar RPA dielectric 
function of the uniform electron gas 

dte-'' J dv(l) (v) tk ■ we-'""-^' (44) 

identifying the excitation spectrum z (k). 

For Z 7^ the modes, z (k, r), depend on r due to the inhomogeneity caused by the ion. 
As an example, consider the solutions with very large z. Expanding e'(r, k; z) to order 1/ z"^ 
gives 

9/, N nir)Vee{k)k'^ ., n{r)Vee{k) ^ n, , 

z^ k, r = ^ ^ ^ ^ + ih\r ^ > ' I3mujl r + .., (45) 

m m 

where = k ■ r/r is the component of k along r. For small k the first term goes to the 
square of the local plasma frequency 

n{r)k^V,,{k) ATin{r)e^ ^ ^ ^^^^ 



and (l45l) simplifies to 



m m 



z^ (k, r) = -ul (r) + lul (r) + .., (47) 

[k\D (r)) 

Here A|, (r) = l/(3mujp{r) is the corresponding local Debye length. Therefore, if k^\ = 
the system supports plasmons with frequencies defined in terms of the local density. This is 
suggestive of a more general "local density approximation"^^ where the dielectric function 
for the uniform electron gas is modified by replacing the uniform density with the actual 
non- uniform density 

e{r,'k;z)^eRPA{ne,k;z) . (48) 

ne=n(r) 

However, this is not correct in general as is evident from fHSl) for k\\ 7^ and the following. 
To continue with the evaluation of fH3l) separate into bound and free contributions 

~ /"°° d [ r 
e{r,k;z) = I + (3n{r)Vee{k) J dte'''- J dv(f) (v) & {v - Vm{r)) e'^-^''^-'^~''^ 

+ ^Mr) I dvQ (vUr) - v) 6 (v?) 6 (u - (r)) e*k-(r(-t)-r) _ ^49) 
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Next, introduce the approximate trajectories of the last section. The analysis is straightfor- 
ward but lengthy so only the final result is given here 

?(r, k;z) = l + pVeeik) {nf{r)If{k, r, z) + n,(r) /^(k, r,^)) . (50) 

The first term of the brackets, proportional to nj(r) = n{r) — ni,{r), is the contribution from 
free states, while the second term proportional to nh{r) is that from bound states (recall 
( !27|) for the bound state contribution n{,(r)). The functions If{k,r,z) and /^(k, r,2;) are 



oo 

oo 







Vm {r)/vo 



If{k, r, z) = kvo dte-^'-^^^ (51) 



Jo 



ik\\Vc{r) sm{uJc{r)t)Jo{kj_Vc (r) sin(a;c('^)^)) 

+ ^Mk±v^ir)sm{uJc{r)t)) . (52) 

The Bessel functions ji(x) and Jo{x) are 

1 1 f^'^ 

jJx) = — (sinx -xcosx) , JJx) = — dcpe'''^""*' . (53) 

X^ 27T Jo 

The free state contribution depends only on the magnitudes k, r while the bound state 
contribution depends on the directions as well, where k\^ and k± are the components of k 
parallel and perpendicular to r. 

Consider the small k, long wavelength limit of this expression for ?(r, k;z). Retaining 
the leading order contributions to fl5T|) and fl52|) gives 



z'^ \ n[r) S^TT \ 



ujlir) ik\\r ribir) 
z^ + u:l{r){k\n{r)f n{r)- ^ ^ 

The limiting forms for these long wavelength excitations as functions of the ion charge 
number Z are 

The local plasmons are recovered for small charge numbers, while for large charge numbers 
the local circular frequencies uj1{r) dominate. The complex coefficient in (155|) implies that 
these excitations are damped. 
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The dielectric function e'(r, k; z) is quite complex, but it is seen that the simple model of 
straight line and circular orbits simplifies this considerably. Further analysis of the collective 
modes for the inhomogeneous electron gas will be given elsewhere. 

VII. DISCUSSION 

A very general description of time correlation functions for electron properties near a 
positive ion has been given by the Markovian approximation (ITSil . The time dependence 
has two contributions, an effective single particle dynamics that dominates for short and 
intermediate times, and a modification of that dynamics due to collective modes. The single 
particle dynamics has a strong dependence on the ion charge number Z, which is due to 
the growing dominance of bound states for large Z. The description is valid for such strong 
coupling conditions, since the Markovian approximation preserves the exact equilibrium 
electron - ion correlations in the effective potential governing the single particle dynamics. 
Similarly the electron - electron potential is renormalized by the exact equilibrium electron 
- electron correlations. 

The description has been illustrated here for the special case of the electric field autocor- 
relation function under the same conditions as have been studied by MD simulations^. As 
the time scales are short, only the effective single particle dynamics has been considered. It 
remains to explore other correlation functions for which the collective modes are expected 
to be more important, and to provide a detailed characterization of those modes via the 
inhomogeneous electron gas dielectric function. The promise for progress in this direction 
is provided by the success of a simple model for the bound and free state dynamics. 

There are several avenues for future directions based on this work: 

1. The analysis here is based on a semi-classical description using a "regularized" electron 
- ion potential. The quality of this type of description can be benchmarked by com- 
parison with a corresponding quantum description of the same effective single particle 
dynamics. 

2. A generalization of the Markovian approximation to a two component electron - ion 
plasma is straightforward. In that case the interest is in the electron dynamics in 
the vicinity of one of the ions. An additional feature is the effect of the dynamics 
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of the ions. This constitutes a kinetic theory for the two particle (electron and ion) 
distribution function. 

3. The spectral line shapes from charged radiators are an important diagnostic tool 
in laser fusion studies. A recent formulation of this problem including all plasma 
charge correlations is expressed in terms of constrained equilibrium time correlation 
functions^. The constraint arises from a specified value of the total ion electric field 
during the dynamical broadening by the electrons. The simple model described here 
provides the potential for practical evaluation of these constrained time correlation 
functions under the demanding conditions of hot, dense matter. The role of charge 
correlations in plasma spectroscopy also has been discussed recently in reference 8. 

4. A corresponding identification of the Markov limit in a fully quantum analysis is 
straightforward but the resulting renormalization of the ion - electron and electron 
- electron interactions by initial correlations is more complicated. Still, the struc- 
ture obtained here of single electron dynamics in the presence of the ion modified by 
collective modes of the dielectric function remains the same. 

5. The attractive distortion of the electron density by the ions, particularly for the bound 
states, is a type of electron confinement. The analysis here can be applied to real traps 
for charged particle confinement (e.g., dusty plasmas near an electrode, ultra-cold plas- 
mas in a laser trap, valence electrons in metallic clusters, electrons in quantum dots). 
Significant differences include complete confinement and relaxing charge neutrality. 
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APPENDIX A: EQUILIBRIUM BBGKY HIERARCHY 

Consider a point ion of charge number Z in an electron gas of average density Ue with 
a positive uniform neutralizing background of density n\, = rif.. The equilibrium structure 
of the electrons in the presence of the ion is given by the one and two particle distribution 
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functions, defined for the equilibrium ensemble by 

fie (no, Vo, Vi) = j dX2. Ax N,pe {^) (Al) 

fiee (rio, r2o, Vo, vi, V2) = Nl j dx3..dxN,Pe (r) (A2) 

where rio = ri — fq is the position of an electron at ri relative to the ion at Tq. This 
position dependence reflects the fluid symmetry (rotational invariance) about the ion. The 
distribution function fie (^"10, vq, vi) obeys the equilibrium BBGKY hierarchy equation 

(vi ■ Vi + Vo ■ Vo + m^^Fie (rio) ■ Vvq + m^^Fei (rio) ■ VvJ fie (no, vo, vi) 
= - j dT2d\2 {niQ^Fie (r2o) ■ Vvq + m'^Fee (r2i) ■ Vvi) fiee (rio, , r20, Vo, Vi, V2) 

+ j dr2 {niQ^Fie (r2o) • Vvq + m^^Fee (r2i) ■ VvJ fie {rio, vo, vi) nf,. (A3) 

The last term on the right side is due to the interaction of the ion and electron with the 
uniform neutralizing background whose density is = ng. Also, Fje (rio) is the force of the 
electron on the ion, F^i = — Fjg is the reaction force of the ion on the electron, and F^e (r2i) 
is the force between electrons. These forces are derived from corresponding potentials 

Fei = -Fie = -^iVei (^o) , F^e (r2i) = -ViKe (rai) . 

The solutions to flA3p have the forms 

fie (ro, Vo, ri, vi) = (f)i{vo)(pe{vi)n (rio) (A4) 

fiee (ro, Vo, ri, vi, r2, V2) = 4>i{vo)(f)e{vi)(f)e{v2)n (rio, r2o) (A5) 
Here 0j(fo) and (f)e{vi) are the Maxwellians for the ion and electron 



^ 3/2 2 



and ?T,(rio), ?^(rio, r2o) are the one and two particle electron number densities (relative to the 
ion) normalized to Ne and N^, respectively. Use of these forms in (lASp gives directly 



= Vo ■ ( - Vin (rio) - /5n (rio) ViVei (rio) 

-P J dr2 [n (rio, r2o) - n (rio) rie] V2Vei (r2o) 



+ vi ■ ( Vin (rio) + (no) ViKi (^o) 

-P j (ir2 (rio,r2o) - (rio)^e] V2V;e (r2i) ^ (A7) 
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Since the velocities of the ion and electron are independent, the following two equations hold 

'n (rio, r2o) 



Vi Inn (rio) = -/^ViKi (^o) - j dr2 
Vi In n (rio) = -/5Vi (no) + f^n^ j 



_ n (rio) ne 
'^(rio,r2o; 



- 1 



- 1 



'^2Vei (r2o) 



V2Vee (r2l) 



(A8) 



(A9) 



n (rio) ne 

This provides two, seemingly independent, equations for the same electron density around 
the ion ('"lo)- Their equivalence imphes 

j dV2 (Fie {r2o) + Fee (?^12)) We (riO, r2o) - U (fio) u^] = (AlO) 

which means that the total external force on the system of two selected particles, the ion and 
the one electron, is zero at equilibrium. In other words, the distortion of riiee (rio, i'2o) is just 
such as to enforce this condition. Both equations (lASP and (]A9P are useful, as illustrated in 
the following two subsections. 

In the remainder of the Appendices and in the text, only the special case of a massive 
ion fixed at the origin is considered. 



1. Screened electric field 



The electric field due to one electron at the ion is defined by 



1 



e (r) = ^VrK. (r) . 



(All) 



The electric field autocorrelation function of ( 123|) depends on the associated statically 
screened field 

1 



es(r) = dr e(r )s (r , r) = e(r) + 



n (r) 



dr'e{r') {n (r, r') — n{r)n (r')) 



(A12) 



This dependence on the two electron density n (r, r') can be eliminated using (1A8I) to get 

1 



e., r 



pZe 



Vi Inn (r) . 



(A13) 



An effective potential is defined in terms of the density n (r) in ([TT 

V.e (r) = -/?-Mnn(r), 



(AM) 



so the screened electric field is given in terms of the gradient of this effective potential 

1 



eslrj = — ViV.e (r) 
Ze 



(A15) 
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2. Hypernetted chain approximation 

A simple and accurate method to determine the density is given by the hypernetted 
chain (HNC) integral equations^. They can be obtained from the usual form for a three 
component plasma of electrons, protons, and ions of charge number Z. Then the limit 
is taken of uniform proton distribution and dilute concentration for the ions of charge Z. 
Instead, it is useful here to state the result as an approximation to flA9|) . Consider the mean 
field limit of no electron - electron correlations 



dX'} 



n (ri,r2 



_n{ri) He 
Then equation flA9|) simplifies to 



V2V;e (^^21) j dY2 



Vi 



n (ri) n (r2) 
n (ri) Ue 



Inn (n) + PV,, (ri) + I3n, j dr, (^^^ - 1^ Ke (rai) 



VaKe (r2i) . (A16) 



(A17) 



In 



n (ri) 



(ri) - Pn, j dr2 - 1^ Ke (rai) • (A18) 

An arbitrary constant has been used to assure the limit n{ri) He when Vei (ri) — > 0. 
Equation ( lAlSI) is an integral form of the Boltzmann - Poisson equation. 

The HNC approximation is similar, but retains electron - electron correlations in the 
absence of the ion 



-(3 dv2 



n(ri,r2 



1 



VaKe (r2i) 



dvo 



n{ri)n{r2) ^ 
n{ri)ne 



V2Cee (r2l) . (Al9) 



n{r\) Tie 

The function Cee {1^21) is the electron direct correlation function defined in terms of the 
electron - electron pair correlation function (without the ion), nee{'r2i), by the Ornstein - 
Zernicke equation 

nee (r) 



Cee (r) = Ke (r) - J dr h^e {t') Cee (|r - T |) , kee {t) 



nt 



(A20) 



The approximation (lAlQp in (1A9P gives the HNC approximation for n (r) 

'^-lVee(r.). 



ln^^ = -/3Kari) + ne / dv2 



(A21) 



This is the same as (lAlSp except that Ke (?"2i) has been replaced by —[3 ^Cee {^21] 
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The electron - electron direct correlation function is determined independently from the 
Ornstein - Zernicke equation flA20|) and the HNC approximation 



In = -f3Vee in) + Ke (t) - c,e (r) . (A22) 



Equations flA20l) - flA22p are the HNC equations used for the numerical calculations presented 
here. 

APPENDIX B: EVALUATION OF C{x,x';t = 0) 

The dynamics of CAsit) is conveniently expressed in terms of the fundamental correlation 
function G{x, x'; t) 

CAB{t) = j dxdx'a{x)G{x,x';t)b{x'), (Bl) 
G{x, x'- 1) = if {x, t) if ix') - if {x'm , f{x) = Y,^{x-x.). (B2) 

a=l 

The initial value G(x,a;'; 0) is easily calculated, with the result 

G{x, x'\ 0) = n (r) [y) [b [x — x') + (f ') n {r') h (r, r')) . (B3) 

In the last equality the two electron correlation function h (r, r') has been identified from 
([7]). Comparison with ([HD shows that h{x) can be written 

n{r)(j){v) J 

This leads to the representation ^ for CAsif) 

CAB{t)= dxn{r)^{v)a{x)b{x,t), (B5) 



with 

b{x,t) = I dx'U{x,x']t)b{x'), (B6) 



U (x, x'; t) = -r^—-T [ dx"G{x, x"- t)G^\x", x'; 0)n{r')^{v'). (B7) 
n{r)(f){v) J 

The inverse of G{x,x'] 0) has been introduced by the definition 

*(.-.',^/<<."G(.,.";0,G-(.'',.';0.) (B8, 
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It is verified that 



G-\x,x';0) 



n (r) (j) (v 



-6 {x — x) — c (r, r') , 



if c (r, r') obeys tlie equation 



c (r, r') = h (r, r') — j dv"h (r, r") n (r") c (r", r') 



(B9) 



(BIO) 



wliicli is a generalization of the Ornstein-Zernicke equation^. 

The formal equation for b{x,t), ([H]), follows from differentiation of ( ]B6I) with respect to 
time and the identification 



C{x, x'\ t) 



dx" {dtU{x, x"; t)) U'^x" , x';t). 



This provides the desired result for identifying the Markovian approximation 

£ {x, x') = C {x, x'; t = 0) = -dtU{x, x'; 0), 



(Bll) 



(B12) 



where the property U {x, x';t = 0) = 6 {x — x') has been used. Equation (1B7P gives finally 

1 



x') 



, , , dx"dtG(x,x";t) G~\x",x';0)n(r'U(v') 
n[r)(j)[v) J t=o 



(B13) 



It only remains to calculate the initial derivative dtG{x,x";t) \t=o to determine C{x,x'). 
To simplify the notation it is useful to denote the force on the electron due to both the ion 
and the uniform positive background by Fq (r) 



dtG(x,x';t) = — dxi6 {xi — x) (vi ■ Vr. + m ^Fo(ri)-Vvi) 
t=o J ^ 



6 (xi — x') n (ri) (t>i) 



+ / dx26 {x2 - x) (j) (f i) {V2) {n (ri, - n (ri) n (r2)) 



j dxidx25 {xi — x)m ^Fee (ri2) ■ Vvi 



X 



{5 {x' - xi) + 5 {x' - 0:2)) {vi) 4> {V2) n (ri, r2) 



+ j dx^5 {x - X3) ivi) {y2) iy^) {n (ri, r2, rs) - n (ri, r2) n (/)) 
(v ■ Vr + m'^Fo (r) ■ Vv) x'; 0) 
-Vv ■ 5{x - x)(t) (v) ttT^ j dY2^ee (|r - r2|) n (r, r2) 



+/?0 {v) {v') V 



Fee(|r-r'|)n(r,r') 



+ j (ir2Fee (|r - r2|) (ri, r2, r') - n (ri, r2) n (r')) 



(B14) 
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The two integrals on the right can be performed using the hierarchy equation (1A9P for n (r) 
and the corresponding next order hierarchy equation for r2(ri,r2). In the current notation 
these are 

j dr^^ee (|r - r2|) n (r, v^) = p-^m (no) - n (no) Fq (r) , (B15) 
J rfraFee (|r - ral) n (n, r^, r') = p-'W.n (n, r')-(Fo (n) + (n - r')) n (n, r) . (B16) 



Then (IA19D becomes 



dtG{x, x\ t) 



t=o 



- (v • Vr + m-^Fo (r) ■ Vv) G{x, x'; 0) 

— (m~^ (/S""^ Vq In n (r)) — m^^Fo (r)) • Vv5 (a;' — x)n (r) (t;) 

{v) (-u') (n (ri, r2) - n (ri) n (r2)) v • 

X Vr In {n (n, ra) - n (n) n (ra)) - Fq (r)) . (B17) 



Next, ehminate the deha function using fIBSP 



dtG{x, x'; t) 



t=o 



- (v ■ Vr + m-^Fo (r) ■ Vv) G{x, x'; 0) 

- (m-i Vr In n (r)) - m^^Fo (r)) ■ VvG(a;, x'; 0) 
+ {m~^ Vr In n (r)) - m"^Fo (r)) 

■Vv</> (w) (f ') (n (r, r') — n{r)n (r')) 
(f ) (f') (n (ri, r2) - n (ri) n (r2)) v • 

X Vr In {n (ri, ra) - n (n) ^ (ra)) - Fq (r)) 

- (v ■ Vr + Vr In n (r)) ■ Vv) G(x, x'; 0) 

+0 (f ) (t;') n (r) n (r) v • V^h (r, r ) . (B18) 



Substitution of this result into ( ]B13[) gives the generator for initial dynamics 
C{x, x') = (v • Vr + m-^ Vr \nn{r)) ■ Vv) 6 {x - x') 



V ■ Vr / dx"(j) {v") n (r") /i (r, r") G'^x" , x'; 0)0 (f ') n (r') 



(v ■ Vr + Vr In n{r)) ■ Vv) 5 ( 



-V ■ Vr 



.(r.r')-/*"„,OMr.r>(r",0 



(v') n (r') . (B19) 



Finally, using (IBlOp the result (fTOl) is obtained 
£(x, x') = (v • Vr + m"^ Vr In n (r)) • Vv) 5 (s - x') - v • VrC (r, r') (i;') n (r') . (B20) 
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APPENDIX C: SOLUTION TO KINETIC EQUATION 

The solution to (fT3|) can be obtained in terms of an effective single electron dynamics by 
direct integration 

b{x, t) = e-^^'bix) - ^ rfre-'^o^'^'V ■ Vr j dr'(3Vee (r, r') /(r, r), (CI) 
where the generator for the effective single particle dynamics is 

£o = V ■ Vr - m- VrHe (r) ■ Vv, (C2) 
and the source term I{r,t) is 

I{Y,t)= f dw(t){v)n{r)h{x,t). (C3) 



The initial condition h{x) is given by ([6]). It is convenient at this point to introduce the 
corresponding Laplace transform 



zti 



b{x,z) = / dte-'%{x,t). (C4) 
Jo 

Then Laplace transformation of (ICip gives the equation z) 

b{x, z) = Gobix) - Gov ■ Vr J dv'(3Vee (r, r') 7(r', z), (C5) 

7(r, z) = j dw(t) {v) n {r)b{x, z), Go = (z + £o)"' (C6) 
An equation for /(r, z) follows from substitution of (ICSP into (]C6I) 

y (iv0 (f ) n (r) Gob{x, ^) ~ J ^'^'^ (^) ^('")^oV ■ Vr J dr'pVee (r, r') -^(r', z) 
7o(r, ^) + y dr'V (r, r"; z) J dr'Vee (r", r') 7(r', (C7) 



where vr fr, r"; 2;) is 



and 



vr (r, r"; z) = -pn{r) j dwcj) {v) Go^ ■ ^rS (r - r") , (C8) 
7o(r, z)=n (r) j d^jcj) {v) Gob{x). (C9) 
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This is an integral equation for J(r, t) which can be written 

j dr'e (r, r'; z) 7(r', z) = 7o(r, z) (CIO) 
The dielectric function e (r, r'; 2;) is defined by 

e (r, r'; z) = S (v - v') - J dv"n (r, r"; z) Vee(r", r'). (Cll) 



(C12) 



With these results ( ICSP becomes 

b{x, z) = Go b{x) - V • Vr y dr"f3Vee (r, r") J dr'e'' (r", r'; z) 7o(r', z) 

The inverse dielectric function is defined by 

dr"e (r, r"; z) e^^ (r", r'; z) = 5 {t — r') = j dr"e^^ (r, r"; z) e (r" , r'; z) . (C13) 

Equation flC12p is the desired solution to the kinetic equation, in terms of the single particle 
dynamics of Qq, since all terms on the right side are now explicit. 

The low frequency limit {z = 0) of e (r, r'; z) has a simple form in terms of the electron 
correlations. First write e (r, r'; z) as 



e (r, r'; z) 



5{r-r') - j dr'n (r, r"; z) Vee(r", r') 
= 6{r-r')+ f3n{r) j dvct> {v) [G^^ - z) Vee(r, r') 

= 5 (r - r') + Pn{r) {v,,{v, r') - z j rfv0 (v) GoVUr, r')^ . (C14) 

where the definition of it (r, r'; z) in flC8l) and Qq in flC6l) have been used. Then taking the 
real and imaginary parts of z going to zero gives 

e (r, r'; z = 0) = 6 {r — r') + /5n(r)Vee(r, r') = 5 (r — r') — n{r)c (r, r') . (C15) 

It follows from the generalized Ornstein-Zernicke equation fIBlOl) that the inverse of 
e (r, r'; z = 0) is 

e-^(r,r';z = 0) = s(r',r), (C16) 
where the static structure factor is defined by 

s (r, r') =5{r-r') + h (r, r') n (r') = 6 {r - r') + h (r', r) n (r) . (C17) 
The high frequency limit {z 00) of e (r, r'; z) also has a simple form 

e (r, r ; z) ^ 5 {r - r ) + ^/5n(r) ^ dv0 (t;) v ■ VrVee(r, r') = 5 (r - r') (C18) 
This implies no screening at asymptotically short times. 
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APPENDIX D: CORRELATION FUNCTIONS 



The Laplace transform of the correlation function ([5]) is 

CAsiz) = / dxn{r)(j) (v) a{x)b{x, z) 



(Dl) 



Substitution of the solution flC12p gives 



CABiz) = / dxn{r)(f) (v) a{x)Qo 



b{x) - V ■ Vr y" dr"pVee (r, r") j dr'e-' (r", r'; z) W, z) 



(D2) 

Considerable simplification occurs for the special case where a{x) = a(r), i.e. it is indepen- 
dent of the velocity, 



Cab{z) = / dxn{r)(f){v)a{r)Qob{x) - 



dTa{r)Pn{r) / dr'" / rfv0 (v) Gqw ■ Vr5 (r - r'") 



X / (ir'Vee (r", r") / dr't ^ {v" ,r' ; z) lair' , z) 



dxn{r)(f) {v) a{r)QQh[x) + j dra{r) J dr'" n {r , r'" ; z) 
X / rfr"Vee (r", r") / dr't'' (r", r'; z) 7o(r', z). (D3) 



where use has been made of the definition ( ]C8I) for tt (r, r'"; z). It follows from ( ICllj) that 
J dr"'n (r, r'"; z) J dv"Vee (r'", r") (r", r'; ^) = e"^ (r, r'; -z) — 5 (r — r') , (D4) 



so (]D3|) becomes 

C^ab(^) = y c/xn(r)0 {v) aiv)gob{x) + ^ c/ra(r) J dr' [e'^ (r, r'; z) - 5 {r - r')] Jo(r', 

(D5) 

dr J dvn (r) (t>) J dr'a{r')e^^ (r, r; z) GoK^) (D6) 

where f lC9P for /o(r', z) has been made explicit, and the dummy labels r, r' have been inter- 
changed. Finally, this can be put in the simple form 



Cab{z) = J drd\n (r) (v) a^(r; z)Qob{^), 



where 



as{r; z) = J dr'a{r')e ^{r',r]z) 
Use has been made of the fact that n (r) (f ) commutes with Cq. 



(D7) 



(D8) 
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1. Dynamic structure factor and field autocorrelation function 

The dynamic structure factor is the autocorrelation function for the electron density, 
corresponding to 

a{x) = 5 (r - q) , b{x) = 5 (r - q') . (D9) 
The correlation function (1D7P in this case is 

C(q, q'; z) = J drdvn (r) (v) e'' (q, r; z) Go (z) s (r, q') , (DIO) 



where s (r, q') is the static structure factor of (]C17|) . 



The electric field autocorrelation follows from flD10|) by integration 



C{z) = J rfqrfq'e(q) ■ C(q, q'; 2)e(q') = j dvdwn (r) {v) e,(r; z)goes{r) (Dll) 

[r; z) = j rfqe(q)e~^ (q, r; z) , e^(r) = j rfqe(q)s (r, ^) = J c?qe(q)e"^ (q, r; 0) 

(D12) 
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